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We report on the observation of cyclotron resonance induced photocurrents, excited by continuous 
wave terahertz radiation, in a 3D topological insulator (TI) based on an 80 nm strained HgTe 
film. The analysis of the photocurrent formation is supported by complimentary measurements of 
magneto-transport and radiation transmission. We demonstrate that the photocurrent is generated 
in the topologically protected surface states. Studying the resonance response in a gated sample 
we examined the behavior of the photocurrent, which enables us to extract the mobility and the 
cyclotron mass as a function of the Fermi energy. For high gate voltages we also detected cyclotron 
resonance (CR) of bulk carriers, with a mass about two times larger than that obtained for the 
surface states. The origin of the CR assisted photocurrent is discussed in terms of asymmetric 
scattering of TI surface carriers in the momentum space. Furthermore, we show that studying the 
photocurrent in gated samples provides a sensitive method to probe the effective masses and the 
mobility of 2D Dirac surface states, when the Fermi level lies in the bulk energy gap or even in the 
conduction band. 

PACS numbers: 


I. INTRODUCTION 

The physics of relativistic Dirac fermions in semicon¬ 
ductors has recently moved into the focus of modern 
research, due to their unique properties and promising 
novel applications, for recent reviews seeA.^. Among 
diverse systems addressed in the literature HgTe-based 
structures represent an extraordinary material class, al¬ 
lowing the fabrication of a high quality material^^— with 
Dirac-like systems of different forms. The latter includes 
topological protected edge and surface states of two- 
and three-dimensional topological insulatorsi^^— , QWs 
with critical thicknes a^^i^^~ — and bulk HgCdTe material 
at the point of semiconductor-to-semimetal transition^. 
Three dimensional TIs based on strained HgTe filmai^ are 
of particular interest. Indeed, in these materials strain 
opens a gap in the otherwise gapless HgTe, which to¬ 
gether with the high quality of the material allows one 
to obtain insulation in the bulk and to study electron 
transport in surface states onlyi^i^i. This differs signifi¬ 
cantly from all other known 3D TI (e.g. Bi 2 Te 3 , Sb 2 Te 3 ), 
where dc electron transport is (almost always) hindered 
by the high residual carrier density in the bulb^Sr— . This 
unique property allows to observe of the quantum Hall 
effecti^ and Shubnikov-de Haas oscillations^! in strained 
HgTe film and, thus to analyse surface state transport in 
TIs. Furthermore, the negligible contribution from bulk 
carriers opens a way to study Dirac fermions in 3D TI by 
cyclotron resonance measured via transmission or Fara¬ 
day effect^Sr.!! and, as we show below, terahertz (THz) 


radiation induced photocurrents. 

Here, we report on the observation of cyclotron res¬ 
onance assisted photocurrents, excited in surface states 
of strained HgTe films. The analysis of the photocur¬ 
rent formation in topologically protected surface states 
is supported by complimentary measurements of the ra¬ 
diation transmission and magneto-transport. We demon¬ 
strate the photocurrent stems from magnetic field in¬ 
duced asymmetric scattering of nonequilibrium surface 
carriers in the momentum space. Our results reveal that 
studying the photocurrent provides a sensitive method 
to probe the effective masses of 2D Dirac surface states. 
Importantly, the method is applicable even for small size 
samples, e.g. gated structures, and for Fermi level ly¬ 
ing in the conduction band, i.e. under conditions where 
transmission experiments are almost impossible. 


II. SAMPLES 

The experiments are carried out on molecular beam 
epitaxy grown high-mobility 80 nm thick HgTe films. The 
HgTe film is sandwiched between thin Cdo. 65 Hgo. 35 Te 
layers acting as capping (top) and buffer (bottom) layers. 
The cross section of the structure is shown in Fig.[TJa). 
The structure is grown on a GaAs substrate with (013) 
surface orientation covered by a thin ZnTe layer and 4 ^m 
thick CdTe layer. The CdTe layer is fully relaxed because 
of its large thickness. The CdHgTe/HgTe/CdHgTe layers 
grown on top adopt the underlying CdTe lattice constant. 
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FIG. 1: (a) Cross section of the investigated structures, (b) 
A X-ray reciprocal space map near the asymmetric (113) re¬ 
flection. The diffracted intensity is plotted on a color coded 
logarithmic scale. The CdTe film is fully relaxed with respect 
to the GaAs substrate, because the line connecting their max¬ 
ima points to the origin of the reciprocal space. The HgTe 
and CdHgTe films are fully strained with respect to the CdTe, 
as their peaks lie on a straight line along the 2 -axis perpen¬ 
dicular to the sample surface, (r. 1. u-reciprocal lattice unit) 




FIG. 2: Experimental setup of (a) the photocurrent and radi¬ 
ation transmission measurements in squared shaped samples 
and (b) photocurrent measurements in cross shaped samples. 


The lattice mismatch between HgTe and CdTe of about 
0.3%, results in an tensile strain in the HgTe film, which 
opens a topological gapi^. The evidence for substantial 
strain in our samples comes from X-ray diffraction mea¬ 
surements, shown in Fig.[Tl(b). Note, that such a sand¬ 
wich design of the structure reduces the influence of dislo¬ 
cations caused by the lattice mismatch^i and allows one 
to obtain high electron mobilities (up to 4x 10^ cm^/V-s) 
together with low residual bulk impurity concentration 
(about 10^®cm“^ as estimated for ungated samples). 

Several kinds of gated and ungated samples have been 
prepared from the same wafer, including ungated square 
shaped samples of 4x5mm size, see Fig.[2Ka), gated (un¬ 
gated) cross shaped structures, see Fig. [2Kb), and gated 
hallbar samples. On the square shaped samples we have 
fabricated eight ohmic contacts in the middle of the 
edges. The sample edges are grown along x || [100] and 
y II [031]. The samples of this kind have been used for 
simultaneous transmission and photocurrent measure- 


FIG. 3: Magneto-transport data obtained in ungated cross 
shaped (a-d), gated cross shaped (e) and gated hallbar (f) 
samples. The panel (a) shows the pxx and (b) the pxy B- 
dependence for T = 4.2 K and T = 40K. (c) Electron Ns and 
hole Ps densities and (d) mobilities {p.e and pp) for different 
temperatures, extracted from the Drude model, (e) Ug de¬ 
pendence of pxx at r = 4.2 K and T = 40 K. The valence and 
conduction band edges are marked by arrows and Ec", 

respectively, (f) Electron and hole densities of bulk 
and and surface states {N^olrac. ^'^^1 N^lx°^). 


merits. In order to study photocurrents and magneto¬ 
transport as a function of the Fermi level position we have 
fabricated gated structures. To avoid insulator leakages 
we have used small area cross shaped structures with four 
ohmic contacts, Fig.|2]|b). The structures of 50 |xm width 
and 1500 pm length have been patterned by means of 
standard photolithography and wet etching. Semitrans¬ 
parent Ti/Au gates of 20 nm/5 nm thickness and 1500 pm 
diameter have been deposited on 100 nm Si 02 and 200 nm 
of Si 3 N 4 dielectric layer grown by chemical vapour de¬ 
position. To monitor the influence of the gates on the 
CR position we have used both, gated as well as un¬ 
gated, cross structures. Gates have also been deposited 
on hallbar samples, which have been used primarily for 
magneto-transport experiments. 

Samples have been characterized by magneto¬ 
transport measurements using standard low-frequency 
lock-in technique with the currents in the range of 0.1 — 
1 pA and with the magnetic field H up to 7 T applied 
normal to the HgTe film plane. The measurements of 
Pxx{B) and pxy{B) for the ungated sample at T = 4.2 
and 40K are presented in the Fig.|3]|a,b). The obtained 
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Pxx{B) shows a large positive magnetoresistance and the 
Hall resistance Pxy{B) exhibits a nonlinear iV-type shape 
around i? = 0 T, which both are typical for a system with 
coexisting electrons and holes. Fitting these traces by us¬ 
ing the classical two-component Drude model we obtain 
electron and hole densities and mobilities'^. The results 
of the fitting are presented in Fig. [21Ic,d). The temper¬ 
ature dependence of electron and hole densities shows 
typical activation behavior at T > 50K and almost sat¬ 
urates below this temperature. The presence of holes at 
low temperatures indicates that the Fermi level is situ¬ 
ated in the valence band while the observed electrons are 
attributed to the surface topological states. Note that 
at high temperatures electrons populate surface as well 
as bulk states. Figure [Ud) shows that both, electron 
and hole, mobilities decrease with increasing tempera¬ 
ture. We emphasize that the electron mobility remains 
as high as 10®cm^/V-s even at T = 40 K. 

Now we turn to magneto-transport data of gated sam¬ 
ples. The Pxx{Ug) dependence obtained at T = 4.2 and 
40 K is presented in Fig. EJe). The low temperature re¬ 
sistivity shows a maximum near Ug ~ —0.8 V. At the 
same gate voltage the Hall resistivity pxy changes its 
sign at B = IT (not shown). This indicates that the 
Fermi level can be tuned from the conduction band to 
the valence band by variation of Ug. Following the proce¬ 
dure described in Ref.H we analysed magneto-transport 
data and identified the gate voltages corresponding to 
the Fermi level lying at the top of the valence band and 
the bottom of the conduction band, see arrows marked 
”if„” and ’’Ac” in Fig. [Sje), respectively. The electron 
and hole 2D density of bulk and surface states as a func¬ 
tion of Ug are shown in Fig.[3Kf). This plot shows that 
B,t Ug ^ OY the Fermi energy is situated in the valence 
band, i.e., holes and surface electrons coexist. At these 
voltages the magnetic field dependence of longitudinal 
and Hall resistance for gated and ungated samples be¬ 
come similar. This indicates that gate fabrication does 
not change the band diagram of the structure qualita¬ 
tively but only introduces a built-in electric field due to, 
e.g., charged impurities and/or defects. The defect den¬ 
sity can be estimated as (1 2) x 10^^ cm^, varying from 

sample to sample. 

For photocurrent excitation and cyclotron resonance 
transmission measurement we apply a cw molecular 
laser— emitting radiation with frequency / = 2.54 THz 
(wavelength A = 118 P-m), / = 1.62 THz (A = 184 pm) and 
/ = 0.69THz (A = 432 pm). The incident power P m 
10 mW is modulated at about 180 Hz by an optical chop¬ 
per. Structures were placed in a temperature variable 
optical cryostat with z-cut crystal quartz windows. The 
radiation at normal incidence is focused onto a spot of 
about 1.5 mm diameter at the center of the sample. The 
spatial beam distribution has an almost Gaussian profile, 
measured by a pyroelectric camera—. The initial linear 
polarization can be transformed into right (cr'*") and left 
(<J~) handed circularly polarized light by applying a A/4- 
plate. A magnetic field A up to 7 T is applied normal 



Magnetic field, B (T) 

FIG. 4: Panel (a): Normalized transmissivity of the square 
shaped HgTe sample at 40 K as a function of magnetic field. 
The data are given for right (cr^, triangles) and left (<j~, 
circles) circularly polarized light. The inset shows the av¬ 
eraged transmissivity traces obtained for different tempera¬ 
tures. Panels (b) and (c) show the photosignal normalized 
by the radiation power {Ux/P) for a'^ and a~ and for linear 
polarized radiation, respectively. Full lines in all panels show 
fits by the Lorentzian function. 


to the film. The photoresponse is measured applying 
standard lock-in technique, either as a photocurrent 
or the corresponding photovoltage Ux,y picked up across 
50 H or 10 MH load resistors, respectively. Parallel to the 
photocurrents we also studied radiation transmission, see 
Fig.HKa). For that a Golay cell detector has been placed 
behind the sample. 


III. EXPERIMENTAL RESULTS 

We begin with the transmissivity data obtained from 
the square shaped sample. The data for r = 40K are 
shown in Fig.UKa) for different polarization states. Excit¬ 
ing the sample with right-handed circularly polarized ra¬ 
diation (cr"^) and sweeping the magnetic field, we observe 
two resonant dips at positive magnetic fields A = 2.6T 
and B = 3.35 T. Upon changing the radiation helicity 
from cr'’' to (t~ the dips appear at negative magnetic 
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FIG. 5: (a) The band structure of (013)-grown 80 nm strained 
HgTe film calculated by k p model. Solid and dashed lines 
show the band dispersion for two perpendicular in-plane crys¬ 
tallographic directions: k || [100] (solid curves) and k jj [031] 
(dashed curves). Calculation are carried out for build-in elec¬ 
tric field Ez = 0.5 kV/cm which results in the energy split¬ 
ting of the top and bottom surface states, shown by thick 
red solid and dashed curves. Inset shows calculated cyclotron 
masses for surface states as a function of energy, (b) Cy¬ 
clotron masses obtained from the transmission measurements. 
Dots and triangles show the results for two CR lines. 


fields. For linearly polarized radiation, the resonances are 
observed for both magnetic field directions (not shown). 
The resonances can be well fitted by a Lorentzian func¬ 
tion. All these features are clear signs that the absorption 
of radiation is caused by cyclotron resonance of electrons. 
The dips in the transmissivity 'T are detected in the tem¬ 
perature range from 4.2 up to 150 K. Below 4.2 K the am¬ 
plitude of both dips is smaller than the noise level of our 
setup. With rising temperature, the two dips increase 
in amplitude and merge together at T above ~ 80 K, 
see inset of Fig. SKa). The temperature dependence is 
due to the increasing carrier density with higher temper¬ 
ature and has been independently checked by magneto¬ 
transport data on the same sample, see Fig. jMc). These 
results are in agreement with findings of Ref. Idol^ where 
double CR structures have been detected in similar sys¬ 
tems by transmission and Faraday effect applying radia¬ 
tion of lower frequencies and attributed to top and bot¬ 
tom surface states. 

From the cyclotron resonance positions given by 


Bcr 


27r/ 


mcR 

e 


( 1 ) 


we determined the correspondent cyclotron masses, 


mcR- At T = 40 K we obtained the masses rric = 0.028 toq 
and me = 0.035mo, for dips at i? = 2.6T and R = 3.35T, 
respectively. To analyse the origin of the resonances we 
carried out k p calculations of the band structure for the 
(013)-oriented 80 nm wide HgTe film sandwiched between 
CdHgTe barriers adapting the procedure described in the 
supplementary materials to Ref.[l^ and Refs. [2i4l The 
details of our calculations for (013)-oriented structures 
taking into account build-in electric fields are given in 
the Appendix. The obtained band structure and effective 
cyclotron masses for various energy levels are shown in 
Fig.lHa) and the inset, respectively. The calculations re¬ 
veal the presence of surface states in the band gap of bulk 
strained HgTe. The Dirac points of the surface states for 
both top and bottom interfaces are deep in the valence 
band which leads to level mixing and deviation of the 
energy dispersion from the linear one, similarly to the re¬ 
sults obtained for (OOl)-oriented structures!^. In (013)- 
grown structures, the energy spectrum is anisotropic in 
the interface plane. To illustrate the anisotropy we plot¬ 
ted the band structure in Fig.[51]a) calculated for two per¬ 
pendicular in-plane directions (solid and dashed lines). 
Moreover, compared to (OOl)-oriented films, we find that 
the spin carried by the surface electrons is tilted out of the 
interface plane, reflecting the low spacial symmetry of the 
(013)-oriented structure, see Appendix. Our calculation 
shows that the masses for surface states are very close to 
the experimentally detected ones and, also in agreement 
with experiment, vary only weakly with temperature, see 
Fig. [5Kb). The detail comparison between the experimen¬ 
tally observed cyclotron masses and the calculated val¬ 
ues require the knowledge of electric field distribution 
and chemical composition profile in real structures and 
is out of scope of the present paper. We note that the 
density of bulk states in the valence band is very high 
due to large effective mass of holes. Therefore, the Fermi 
level in ungated samples and in gated samples at neg¬ 
ative voltages is efficiently pinned to the valence band 
top. Thus, we conclude that the resonances are caused 
by the top/bottom surface states. This conclusion is also 
supported by magneto-transport data indicating a van¬ 
ishing amount of bulk electrons. Moreover, according to 
the calculation, cyclotron resonance of the bulk carriers 
is expected at substantially larger magnetic fields corre¬ 
sponding to a mass of about 0.07mo. 

Next we demonstrate that the CR, so far observed in 
transmission, results in a photocurrent, whose magnetic 
field and polarization dependencies fully reproduce the 
peaks positions and shapes of the CR absorption. The 
magnetic field dependencies of the corresponding photo¬ 
signal are shown for circularly or linearly polarized radi¬ 
ation in Fig.SKb) and (c), respectively. The coincidence 
of the transmission and photocurrent resonance positions 
indicates that the current is also excited in the surface 
states. For a more detailed study of the photocurrent we 
used cross shaped samples, which are too small to be used 
for transmission measurements, but allow fabrication of 
semitransparent gates. Figure [5] shows the resonant pho- 
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FIG. 6: Normalized CR radiation induced photocurrent mea¬ 
sured for different radiation frequencies in ungated cross 
shaped samples. Note that the data are j/-shifted by 0.2 for 
and 0.4 for / = 1.62 and 0.69 THz, respectively. Full lines 
show fits by the Lorentzian function. Inset shows the first, 
Bern (squares), and second, Bcr 2 (circles), resonance posi¬ 
tions as a function of frequency. Full and open triangles show 
the data of Ref.|^ for top and bottom surface states measured 
in similar strained HgTe film of 70 nm width. 


tocurrent obtained for different radiation frequencies in 
an ungated cross shaped sample. Comparison of the data 
for / = 2.54 THz obtained for square and crossed shaped 
samples, shows that magnetic field positions of the res¬ 
onances stay almost the same. Notice that the relative 
strength of the first, Bern, and second, Bcr 2 ^ resonances 
substantially change with decreasing frequency. While 
at / = 2.54 THz the second resonance is substantially 
stronger than the first one, for / = 0.69 THz the sit¬ 
uation reverses. The inset in Fig.jS] demonstrates that 
the position of both resonances linearly scales with the 
frequency. The frequency dependence of the CR is in 
line with transmissivity data of Ref.H^ of top and bot¬ 
tom surface states, obtained for a strained HgTe films of 
70 nm width, see triangles in the inset in Fig. [5] In order 
to provide additional support for the conclusion that the 
observed resonant photocurrent and CR stem from two- 
dimensional surface states, we carried out measurements 
with the magnetic field tilted by an angle 0. Figure[7] 
shows that the current is caused by the magnetic field 
component B^, normal to the surface. Indeed, calculat¬ 
ing from the peak position Bqr in Fig. [71 the value of 
Bf^ = Bcr * cos(0), we obtain that B^^ does not de- 



FIG. 7: (a) Magnetic field dependencies of the photocurrent 
measured for magnetic field tilted by an angle 0. The inset 
shows the resonance field, B^^ = Bcr*cos{Q), as a function 
of angle 0. 



Magnetic field, B (T) Magnetic field, B (T) 


FIG. 8: (a) Magnetic field dependence of the photoenrrent 
measured for gated cross sample and different gate volt¬ 
ages, Ug. Fnll lines show fits by the Lorentzian function. Note 
that the data for various Ug are y-shifted for clearness. The 
shifts are: -0.1 V for Ug — -0.5 V; OV for f/g = OV; 0.15 V 
for Ug = 0.5V; 0.35V for Ug = 1.5V; 0.9V for Ug = 3V; 
3.5V for Ug = dV. (b) Gate voltage dependencies of effec¬ 
tive masses determined from the position of CR resonances 
at magnetic field strength B > 2.5 T. (c) Magnetic field de¬ 
pendence of the photocurrent measured at low temperature 
and high magnetic fields. 

pend on the angle 0, see inset. Moreover, for an in-plane 
magnetic field the photocurrent vanishes (not shown), in¬ 
dicating the two-dimensional nature of the carriers^. 

By using gated cross shaped samples we study the pho¬ 
tocurrent and CR as a function of the Fermi energy, see 
Fig.Ha). The data show that at negative bias voltages, 
for which the Fermi energy lies in the valance band, the 
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Magnetic field, (T) 

FIG. 9: Magnetic field dependence of the photosignal mea¬ 
sured for two in-plane directions (x and y) obtained for gated 
cross sample #2 {Ug = 0). 


photocurrent resonances vanish. For zero gate voltage 
a resonance at B = 3.3 T is clearly resolved. With in¬ 
creasing gate voltage we additionally detected a reso¬ 
nance at smaller magnetic field, see data for Ug = 0.5 V 
in Fig.lH^a). The positions of the resonances are simi¬ 
lar to the ones obtained for ungated samples, see Figs.|3] 
and m and correspond to the response of the top and 
the bottom surface states of the film. At f/^ > 1V we 
detect an additional resonance at B Ki 2.4 T. Note that 
for Ug > lY and T = 40K the Fermi energy lies close to 
the bottom of the conduction band resulting in a notice¬ 
able electron density in the bulk, see Fig.[3lje). Moreover, 
the high gate voltage applied to the film may form an in¬ 
version layer for bulk electrons changing their dispersion. 
The origin of this additional resonance is out of scope of 
the present paper. At even higher gate voltages Ug > 2Y 
and T = 4.2 K we detected another resonance at much 
stronger magnetic fields B Ri 6.3 T, see Fig.[5Jc). The po¬ 
sition of this resonance does not depend on the gate volt¬ 
age and the corresponding effective mass is m = 0.07mo. 
This mass matches the one calculated for bulk electrons, 
so that the resonance can be attributed to the CR of 
bulk electrons. We note that at T = 4.2 K the surface 
state resonances are superimposed by Shubnikov-de Haas 
oscillations, which makes the CR analysis almost impos¬ 
sible. 

From the resonance positions of the surface state re¬ 
sponse we determine the gate voltage dependence of the 
effective masses, see Fig.jSJb). While the mass corre¬ 
sponding to Bqri remains almost independent of the 
gate voltage, the Bcr 2 slightly rises with increasing Ug. 
The resonance at B^ ~ 3.3 T {mcR ~ 0.036mo) is at- 
tributet to the bottom surface, whose electron density 
is not much influenced by the gate voltage. In fact it 
is known from the magneto-transport measurements, see 
Fig-El^f); that the filling rates dNs/dUg for top and bot¬ 
tom surface states are considerable different because of 
electrostatic screening of the bottom surface by the top 


one. Moreover the electron density at the top surface 
rises with the gate voltage increase, see Fig.[31Jf), which 
leads to the increase of the cyclotron mass. By contrast 
the density of the bottom electrones is almost indepen¬ 
dent on the gate voltages, see Fig.[3jf). From the full 
widths at half maximum of the CR measured in trans¬ 
mission in square shaped ungated samples, see Fig.|Tl[a), 
we also estimate the mobility of the surface electrons, be¬ 
ing /i = 5.5 • 10^ cm^/V-s. This value fits well to the mo¬ 
bility obtained from magneto-transport in ungated cross 
shaped samples, see Fig.[31[d). Finally, we discuss the 
photocurrent behavior when reversing the magnetic field 
direction. The magnetic field dependence of the pho¬ 
tocurrent measured for two in-plane directions {x and y) 
is shown in Fig.[S] First of all this plot indicates that 
both resonances behave independently of each other. In¬ 
deed, while the resonance for the photosignal Ux at Bcri 
is even in magnetic field, the one at Bcr2 is odd. This 
result agrees well with the above conclusion, that the 
photocurrents are generated by two independent surface 
states (top and bottom). For the other direction this sit¬ 
uation just reverses: Now the resonance at Bcri is odd 
and Bcr2 is even in magnetic field. 

This H-dependence of the photocurrent can be under¬ 
stood from the phenomenological description based on 
symmetry arguments. The appearance of magnetic field 
induced photocurrents, excited by homogeneous radia¬ 
tion, is well known for two dimensional systems lack¬ 
ing spatial inversion symmetry, for review see Ref.H^H^. 
For free carrier absorption, being accompanied by elec¬ 
tron scattering from phonons or static defects, the cur¬ 
rent is caused by the asymmetry of electron scattering in 
momentum space induced by an external magnetic field, 
when the electron systems is driven out of thermal equi¬ 
librium by light absorption^. At cyclotron resonance, 
the electron-photon interaction is enhanced, which leads 
to a resonant behavior of the photocurrent. 

Considering the generation of asymmetric electron dis¬ 
tribution in the momentum space by cw radiation, the 
distribution evolution in the magnetic field and its decay 
due to scattering processes, one can write the Boltzmann 
equation for the distribution function fp 

e[v X B\ ■ — gp ^ , (2) 

where p and v are the electron momentum and velocity, 
respectively, gp is the generation rate of electrons in the 
state with the momentum p, t is the momentum relax¬ 
ation time, and {fp) is the distrubution function averaged 
over the momentum direction. Multiplying Eq. m by ev 
and summing up over p we obtain the equation for the 
photocurrent density j = J2p ^'^fp in the absence of an 
in-plane bias 

j X Wc = G - - , (3) 

r 

where u)c = oJc = eB^/mcR is the cyclotron fre¬ 
quency, z is the unit vector pointing along the z axis. 
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and G = ^'^9p is the rate of current generation. 

The direction of the vector G is determined by the 
point group symmetry of the structure, radiation polar¬ 
ization state and magnetic field orientation. In topolog¬ 
ical states on (013)-oriented surfaces with no nontrivial 
symmetry elements (Ci point group), the current gener¬ 
ation becomes possible even at normal incidence of ra¬ 
diation and magnetic field perpendicular to the surface. 
In this particular geometry the generation rate G can be 
presented in the form 

G = , (4) 

where / is the radiation intensity, r](u!) is absorbance 
and 7 is the vector determined by the magnetic field in¬ 
duced asymmetry of electron scattering processes. For Ci 
point group the direction of the vector 7 is not forced to 
a certain crystallographic axis by symmetry arguments 
and may depend on temperature, gate voltage, radia¬ 
tion frequency, etc. Such a specific anisotropy has been 
demonstrated for photogalvanic currents excited by ter¬ 
ahertz radiation in (013) HgTe quantum wells at zero 
magnetic field^. For two perpendicular in-plane direc¬ 
tions the solution of Eq. ([3]) has the form 


jx — 

jy — 


lx - UJcT -fy 

I -I- (wcr)^ 

Jy + UJcT Jx 
1 + (WcT)^ 


■Ir]{uj)B^ . 


(5) 


Equation ([5]) shows that the current (i) has a resonant 
behavior a.t lu = uic caused by CR in the absorbance 
r]{uj) and (ii) contains contributions even and odd in the 
magnetic field, both in agreement with experiment.The 
formation of photocurrent for positive and negative mag¬ 
netic fields Bz is illustrated in Eig. (TUI sketches the gen¬ 
eration rate vector G for positive magnetic fields. The 
resulting photocurrent j deviates from G due to the Hall 
effect. For opposite magnetic field the vector G and the 
Hall rotation axis reverse as shown in Fig. fTUI b'). Conse¬ 
quently, for X and y projections we obtain both odd and 
even components of the total photocurrent, see Fig. [9] 


IV. SUMMARY 

To summarize, our results demonstrate that excitation 
of 80 nm thick strained HgTe films by terahertz radiation 
results in cyclotron resonance induced photocurrents gen¬ 
erated in the topologically protected surface states. The 
effect emerges due to magnetic field induced asymmetric 
scattering of nonequilibrium surface carriers in the mo¬ 
mentum space. Our study, accompanied by complimen¬ 
tary transmission and magneto-transport measurements 
show that the observed CR assisted photocurrent pro¬ 
vides a sensitive method to probe the effective masses 
of 2D Dirac surface states. Importantly, the method is 
applicable even for small size samples, e.g. gated struc¬ 
tures, and for the Fermi energy lying in the conduction 



FIG. 10: The vectors of photocurrent generation rate G and 
the steady-state photocurrents j for two opposite magnetic 
field directions, panels (a) and (b). The direction of photocur¬ 
rent j in the magnetic field is declined from the direction of 
G by the Hall angle arctan(tJc'r). Figure illustrates that, for 
the given parameters, jx changes its sign while jy converses 
its sign when the magnetic field direction is reversed. 


band, i.e. conditions where transmission experiments are 
almost impossible. 


V. APPENDIX 


We calculate electron spectrum in low-symmetry 
strained CdHgTe/HgTe/CdHgTe structure on CdTe in 
the framework of 6 -band k-p model which is relevant for 
narrow gap semiconductors. We consider the class of 
(OZh)-oriented films (Z and h are integer numbers) which 
includes (Oil)-, (012)-, and (0I3)-oriented structures. In 
the basis of Bloch amplitudes of the Fe and Fg bands, 
the wave function has the form 


'S{p,z) 


/ '0r6.-ti/2 \ 
V'r6.-i/2 
V'rs.-1-3/2 
V'rs.-1-1/2 
V'r8,-i/2 

\ V'r8,-3/2 / 


exp(/fc|| • p). 


( 6 ) 


where p = {x,y) is the in-plane coordinate, z is the 
growth axis, tpj are the envelopes, fc|| = {kx,ky) is the 
in-plane wave vector. 

The effective 6 -band Hamiltonian is given by 


H = 


CV 



(7) 


where the blocks Hcc and describe the conduction 
(Fg) and valence (Fg) intraband contributions, and the 
block Hcv describes the band mixing. The block Hcc is 
given by 


Hcc — I2x2 


Ec{z) 


T? k[2F{z) + l]fc 
2mo 



( 8 ) 


where 12 x 2 is the 2 x 2 identity matrix, Ec{z) is the con¬ 
duction band profile, k = {k^, ky, —ididz)-, F{z) is a pa¬ 
rameter accounting for contribution from remote bands. 






















Sc is the Fg-band deformation potential constant, and e 
is the train tensor. The block has the form 


Hr.ii — 


^ -^Pk+ 
72 


PK 


76 


Pk- 


0 


V 


0 


- ^Pk^ 

76 


2 1 
3^^^ 71 


+ 


(9) 

where P is the Kane matrix element and k± = kx iky. 
The block Hyy is given by 


Hyy = Ey (v) T T T T iF 17 , (10) 


where Ey{r) is the valence band profile, 

(a) 

and Etgp are the isotropic and anisotropic parts of the 
Luttinger and Bir-Pikus Hamiltonians, 



2too 



5 

7l + 2^2 


k + 2(J • k)'-f2{J ■ k) , 
( 11 ) 


^BP = Tre- 6 ega - ^ {Fg, Jp}seap , 

^ ^ a a^/3 

( 12 ) 

7 i, 72 and 73 are contributions to the Luttinger param¬ 
eters from remote bands, J is the vector composed of 


J 


8 

the matrices of the angular momentum 3/2, {Ja, Jf 3 }s = 
(Fg J /3 -f J /3 Ja)/2, a, b and d are the Tg-band deformation 

fa) 

potential constants. The explicit form of the terms iF/ ' 

fa) 

and Hgp related to the cubic anisotropy of host crystals 
depends on the coordinate frame used. 

In the coordinate frame relevant to (0^/i)-oriented 
structures: x || [ 100 ], y || [0F7], and z || [ 0 /F], the terms 
FF 7 ^ and FF^p assume the form 



7 (eV) 

7i 72 

73 

2mo{P/nf (eV) 

HgTe 

CdTe 

-0.303i^) 

1.606C) 

4.1® 0.5^^^ 1.3^^) 

1.47® -0.28® 0.03® 

18.8(^> 

18.8® 


F 

a (eV) b (eV) d (eV) 

H. (eV) 

HgTe 

CdTe 

0 ( 1 ) 

-0.09(^> 

-0.13'^) -1.5® -8.0® 

0.756^®^ -1.0® -4.4® 

-3.82® 

-2.687® 


ao (A) 

cii (Mbar) ci 2 (Mbar) 

C 44 (Mbar) 

HgTe 

CdTe 

6.46<®^ 

6.48® 

0.597® 

0.562® 

0.415^^^ 

0.394^^^ 

0.226^^' 

0.206® 


TABLE I: Parameters of HgTe and CdTe. i^^Ref. Iisl . 

(^^Ref.y, i^^Ref.S WRef.il, (^'Ref.lH, (^’Ref.lH. 


^(a) 


2r 

mo 


f 

\ 1 

A xJy\s{ri^ ^2)^x^y {'Jx^z\ 0^2; ^z\s ^x 

{JyJz}s COS 20 H- ^—- sin 20 


(13) 


^(a) _ 
ripp — 


-2 


73 


-b 


X {(73 - 72 ), 7 }sfcyCOs 20 

{J.Jy}sexy + { Jx + 


7(73 - 72)7 - (73 - l 2 ) k . 


- .P 

{JyJz}sCOs29 -\—^ sin 20 


sin 20 


eyz cos 20 


- e 


— sin 20 


where 0 = arctan( 7 F) is angle between the growth direc¬ 
tion [Olh] and the [001] axis, 0 ~ 18.4° for (013)-oriented 
structures. 

HgTe and Cdo. 65 Hgo. 35 Te layers adopt the in-plane lat¬ 
tice structure of CdTe buffer which leads to the in-plane 
strain of the layers described by the strain tensor com¬ 


r 


ponents = Syy = ocdTe/oo - 1 , where acdTe is the 
lattice constant of CdTe and uq is the equilibrium lat¬ 
tice constant of the considered layer. Out-of-plane com¬ 
ponents of the strain tensor in each layer can be found 
by minimizing the elastic energy. Such calculations for 
(OZF)-oriented structures yield 


cfi + 2 cii(ci 2 - C44) -I- ci 2 (- 3 ci 2 + IOC 44 ) - (cii -I- 3 ci 2 )(cii - C12 - 2C44) cos 40 
-cfi - 6C11C44 -I- Ci 2 (ci 2 + 2C44) -I- (cii -I- Ci2)(cii - C12 - 2C44) cos 40 


(cii -I- 2 ci 2 )(cii — C 12 — 2 C 44 ) sin40 




— cfi — 6 C 11 C 44 + Cl 2 (ci 2 -I- 2c 44) + (cil -I- C 12 )(cii — Cl 2 — 2c44) COS 40 


(14) 
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^xz = 0, where cn, C12, and C44 are elastic constants. 

We have performed calculation of the electron spec¬ 
trum and wave functions for the parameters listed in Ta¬ 
ble I. The Schodinger equation with the matrix Hamil¬ 
tonian 0 is solved numerically by expanding the en¬ 
velope functions '4^j{z) in series of the harmonic oscilla¬ 
tor functions following the procedure described for (001)- 
oriented structures in Ref.l^. The calculated spectrum 
for the (013)-oriented structure is shown in Fig. [5] The 
wave functions of states emerging in the band gap of bulk 
strained HgTe are shown in Fig. [TTJ The states are lo¬ 
calized at the film interfaces and predominantly formed 
from the heavy-hole states. Despite the fact that 6-band 
k p model used in the calculations neglects bulk inver¬ 
sion asymmetry of HgTe and the corresponding Dressel- 
haus spin-orbit interaction, we have found that the spin 
carried by surface states in (013)-oriented structures for 
kx ^ 0 has an out-of-plane component. The out-of-plane 
polarization defined by 

p _ /(IV'rg,+3/2p - |#g,-3/2n dz 

/(IV'r8,-h3/2P + |'0r8,-3/2P) dz 

is about 0.1 for the surface states presented in Fig. [11] 
The emergence of the out-of-plane spin component is 
attributed to the cubic anisotropy incorporated in the 
Hamiltonian ([Tj) together with inversion symmetry break¬ 
ing at surfaces. Indeed, the 6-band k p Hamiltonian of 
bulk crystal corresponds to the Oh point group contain¬ 
ing the spatial inversion. At a surface, the symmetry 
z — z is broken which leads, for (013)-oriented surface, 
to the reduction of system symmetry to the Cg point 
group containing the only non-trivial element x ^ —x. 
Such a symmetry reduction enables the emergence of out- 
of-plane spin polarization for surface states with ^ 0. 
We note that Dresselhaus spin-orbit interaction steming 


from bulk inversion asymmetry of the host zinc-blende 
structure, which in not included in the model Hamilto¬ 
nian, will lead to an out-of-plane spin polarization for 



FIG. 11: Wave functions of edge states in (013)-grown 80 nm 
strained HgTe film. Calculations are carried out for the elec¬ 
tron energy £1=11 meV, the wave vector components > 0 
and ky = 0, and the build-in electric field = 2 kV/cm. 

surface states with ky ^ 0. 
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